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We study gravitational creation of light fermions in the presence of classical scalar metric pertur- 
bations about a flat Friedmann-Lemaitre- Robertson- Walker (FLRW) background. These perturba- 
tions can be large during preheating, breaking the conformal flatness of the background spacetime. 
We compute numerically the total number of particles generated by the modes of the metric per- 
turbations which have grown sufficiently to become classical. In the absence of inhomogeneities 
massless fermions are not gravitationally produced, and then this effect may be relevant for abun- 
£\j , dance estimates of light gravitational relics. 
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I. INTRODUCTION 

(N 

\Q \ It is well known that particles are created in an expanding Universe . Pioneering work by Parker (^] highlighted the 
• creation of nonconformally-invariant particles even in flat Friedmann-Lemaitre- Robertson- Walker (FLRW) Universes. 
The extension to anisotropic Bianchi cosmologies by Zel'dovich and Starobinsky ||, showed that massless, conformally 
coupled, scalar particles are created due to the breaking of conformal invariance which, in four dimensions, is signaled 
by a non- vanishing Weyl tensor. In the presence of small, inhomogeneous metric perturbations, Horowitz and Wald Q| 
obtained the vacuum expectation value of the stress tensor for a conformally invariant field and showed that particle 
production only occurs if (T^ v (x)) is non-local (i.e., has non-zero contributions with support on the past-null cone of 
\ the event x). 

The number of particles created by inhomogeneities can be computed by a perturbative evaluation of the S- 
Q H ' matrix . In the inhomogeneous case, the occupation number of produced particles is composed by three parts, 
viz, the zeroth-order contribution due to the homogeneous expansion, a first-order (in the perturbation amplitude) 
1 part arising from the interference between 0- and 2- particle states, and a second-order contribution which comes 
, from the interaction between nonzero particle states. In the case of massive or nonconformally coupled scalar fields, 
'kJi ' ^ ne last two terms (respectively linear and quadratic in the metric perturbations) typically give a contribution small 
rS | compared with the first one. However, in the massless and conformally coupled cases, the first two terms vanish in 
rigid, exactly FLRW backgrounds, and only inhomogeneity contributes to the gravitational particle production. 

Such theoretical studies were historically of interest both because they focus on the interplay between quantum 
field theory and gravity and because of their relevance to Misner's "chaotic cosmology" program || in which arbi- 
trary initial inhomogeneity and anisotropy were to be damped to acceptable levels due to particle production. The 
non-renormalisability of standard quantum gravity, the successes of string theory and the dominance of cosmological 
inflation have since removed much of the motivation for studies of particle production from inhomogeneity, though 
research on non-equilibrium issues, generalized fluctuation-dissipation relations and effective Einstein-Langevin equa- 
tions still continues ||. 

Inflationary cosmology seems to be less affected by the complex aspects of gravity beyond the semi-classical realm. 
For inflation to start, a homogeneous patch larger than the Hubble radius, H , is needed which may require some 
fine-tuning, especially for inflation at low energy scales. However, once inflation has begun, the no-hair conjecture 
fiof should ensure that pre-existing inhomogeneities are driven to zero. The quantum metric perturbations generated 
during inflation can then be tuned to have small amplitudes ~ 10 -5 and unless one wants to compute 4-point correla- 
tion functions which include graviton loops, comparison with the anisotropies in the Cosmic Microwave Background 
(CMB) is straight-forward. 

However, nonperturbative production during the coherent oscillation of the inflaton fields pd|-p^[, dubbed pre- 
heating Jl2j, can alter the simple picture of small amplitude metric perturbations. Preheating therefore provides a 
particularly useful arena for examining some aspects of the boundary between semi-classical and quantum gravity. 
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Preheatin g ca n lead to the growth of metric perturbations, which in turns is known to stimulate scalar particle 
creation |16H18[| and the production of seed magnetic fields Jl9|-|23|. The amplification of metric perturbations may 
also yield runaway instabilities due to the negative specific heat of gravity, which leads to an interesting possibility of 
primordial black hole (PBH) formation J2j,^| . It can also amplify metric perturbations on very large scales |]l9| , ^6| -p0| . 
The growth of metric perturbations, in the case where a PBH does not form, is limited by backreaction effects which 
are explicitly nonlinear. 

Hence on certain scales the Weyl tensor can actually be quite large. This means that there can be significant particle 
production even if the field is governed by a conformally invariant equation of motion, the classic examples being 
photons, massless fermions and conformally coupled scalars. In particular, if the effective mass of gravitational relics 
such as modulini and gravitini is small enough during reheating, then their main source of gravitational production 
will come from metric perturbations. Since stringent upper bounds hold for the abundances of these species, it is 
important to verify that creation from inhomogeneities (if sufficiently amplified at preheating) does not overcome this 
threshold. Our analysis suggests that this does not occur in the simplest models of chaotic inflation, at least when 
rescattering is neglected. 

Particle creation via inhomogeneities may also be relevant for alternative scenarios such as the pre-big-bang model 
pf| , where homogeneous gravitational production of moduli and gravitinos is known to be an issue 
Universe goes through a high-curvature phase where higher-order a' corrections are important. 

The paper is organized as follows. In section II we review the formalism for computing the number of produced 
fermions, while in section III we discuss the evolution of metric perturbations during preheating in some specific 
potentials. Sections IV and V present our numerical results and conclusions respectively. 



II. FORMALISM FOR PRODUCTION OF FERMIONS 



In this section we review the formalism for the production of (Dirac) fermionic quanta by cosmological inhomo- 
geneities, summarizing the calculations of [^-[7). We consider perturbations around the a FLRW background 



ff/iK = a 2 (v) ( r l^u + h^ v ) 



(2.1) 



where a(r\) is the scale factor and rj the conformal time. We will also use proper time, f, related to conformal time by 
drj = dt/a and choose to consider only scalar metric perturbations in the longitudinal gau ge J33| , |34[ []. 

At linear order minimally coupled scalar fields do not induce an anisotropic stress p7j and hence the metric 
perturbations are characterized by a single potential: /i M „ = 2$<5 M „ , viz: 



ds 2 = a 2 (l + 2$)<V - a 2 (I - 2$)8 i jdx i dx> . 
The action of a fermionic field ip with mass m in the background (|2.2|) is, to first order in $, 
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S 



d 4 xa 3 (1-2$) V ^ 



3 da 



3d$ 



(1-*) 7°^ + --t- (1-$) 7°-7T-^7° 



2a drj 



2 dr] 



(1 + $) fdi--Yd^ 



d x tp 



~ am j tp 

(i-f°d v +iYd l -am) -3$ 7 °ft ? -$7 i ft 



3d$ 
~2d7\ 



y ft* 



d A x [C +d ($)] 



(2.3) 



where we have defined ip = a 6 1 2 ip and Cq, Ci denote the homogeneous and interaction Lagrangians 



= a^ 2 ', 

The number density of produced particles is easily computed in the interaction picture, where the evolution of the 
operators is just determined by the homogeneous expansion of the Universe, and the states evolve according to the 
small inhomogeneities in the metric. The field ip can be decomposed in the standard way 



1 Vector perturbations will induce magnetic fields via the induction equation [|3^] but we do not consider them since inflation 
drives vector perturbations to zero. We also neglect tensor (gravitational wave) modes since they are not strongly amplified at 
preheating. 
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T^T^* br (k, V) a r (k) + v r (fc, r,) b\ (-k)] , (2.4) 

where, as usual, the anti-commutation relations are 

{a r (k), at(k')} = {b r (k), bl(k')} = 5™ (k - k') S rs . (2.5) 

The Fock space is built at the initial time rji starting from the vacuum state defined by 

a r (fc) |0) = b r (fc) |0) = . (2.6) 

For massless fermions (the conformally coupled case) , the expansion of the Universe factors out of the free action 
/ d A x Co and the spinors u r and v r evolve as in Minkowski spacetime. As a consequence, the operators a r W , 6 r (t) are 
the physical annihilation (creation) operators at all times. Introducing a mass term breaks the conformal symmetry, 
and the physical annihilation/creation operators are defined through the Bogolyubov transformations (see |3(| for 
details) 

a r (fc, rf) = a r {k, rj) a r (k) — f3 r * (fc, rj) b r ' (— k) , 

SJ (k, rj) = (3 r (k, rj) a r (k) + a r * (k, r,) bj (- k) . (2.7) 
The two Bogolyubov coefficients have initial conditions a (r/i) = 1 , (3 (iji) — and evolve according to J|7j 

a r = — — ^- e J r , (3 r = e J 'a r , (2.8) 



with u> = yW + aPrnr . These equations preserve the normalization \a r \ 2 + \/3 r \ 2 — 1 , which holds in the fermionic 
case. 

In the interaction picture, the evolution of the initial vacuum (zero particle) state |0) is determined by the interaction 
Lagrangian Cj(Q) . At first order in $ we have 

|V) = |0) + i / d 3 k d 3 k' \k, r ; k', s) (fc, r ; k', s\S\0) , (2.9) 



2 . 

with the 5-matrix element (T stands for time ordering) 

(k,r; fc'^IS'lO) = iT(k,r; k',s\d ($) |0) . (2.10) 

Rigorously speaking, the occupation number can be computed only in the asymptotic future, 77 y , with vanishing 
perturbations &(r)f) = 0. Identical results are obtained for particles and antiparticles, so we concentrate on the 
former. The expectation value of the number operator N = (2 tt a) 3 J d 3 pd\. (p, rjf) d r (p, r]f) in the state \ip) is given 
by the sum of three terms, No + N± + N2 , which are of zeroth, first, and second order in $ , respectively [§-0, 

^o-ttt^tt [ d 3 k(0\al(k)a r (k)\0) = [ d 3 k\f3 k \ 2 , (2.11) 



(27ra) 3 J x 1 ry ' ' y ' 1 ' (27ra) 3 

(2™^ 
1 

4 (2na 



N x = - - 1 Nq [ d 3 k Re [a r (fc) (3 r (fc) (k, r ; - fc, r\S\0)} , (2.12) 
^2= ,^ , 3 [ d 3 kd 3 k'\(0\S\k,r; k',s)\ 2 [| Qr (fc)| 2 + |&(-fc')| 2 + l] . (2.13) 



The zeroth-order term (2.11) is the well known expression arising from the homogeneous expansion of the Universe. 
V denotes th e vo lume of the Universe at late times, when the perturbations can be neglected. The first-order 
contribution ( 2.1 2j ) comes from the combined effect of expansion and inhomogeneities. These two terms vanish for 
massless fermions, which, as we remarked, are conformally coupled to the FLRW background (more explicitly, we 
notice that f3 r (fc, 77) = in this case, as ca n be seen from eqs. (2.8)). 

In the massless case, only the last term (2.13) contributes to the particle production. Furthermore, in this limit N2 
acquires the particularly simple form || 
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in terms of the Fourier transform 



(jabcd ^ = d i xe ipx C abcd (x) 



(2.15) 



of the Weyl tensor C abcd . 

For the metric (2.2), one finds (see also t| f° r useful intermediate steps) 
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15 (Ztt) a J Jo •'IpKpo 



drie ipor >$ (p,r?) 



(2.16) 



where $ is defined by 



$ (rj , x) 



(27T) 
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*(p, »7) e ip ' : 



(2.17) 



In these expressions, $ should be regarded as a purely classical perturbation. The normalization as well as the 
classicality condition for <E> will be discussed in the subsequent sections. 

III. EVOLUTION OF SCALAR METRIC PERTURBATIONS 

A. Linearized equations and analytic solutions for metric perturbations 

We now discuss the evolution of scalar perturbations $ during inflation and reheating in simple models of chaotic 
inflation |38|, denoting the minimally coupled inflaton field by 0, with potential 



V = -m 2 
2 



or 



(3.1) 



In subsection 1IIB 5| we will then discuss a model where two fields are present. 

For the potentials (3J), the COBE normalization of the Cosmic Microwave Background Radiation |3!| requires 
the coupling constants be m ~ 10 _6 M p for a massive inflaton and A ~ 1CP 13 for the quartic case. Decomposing the 
inflaton as cf>(t,x) — * <f)(t) + 5<f)(t, x), the background equations are given by 



H 2 



8tt 
3Af p 2 



V 



3H<f) 



dV 







(3.2) 
(3.3) 



where dots denote derivatives with respect to physical time. In the numerical simulations which we present in the next 
section, we will implement the backreaction of field fluctuations within the Hartree approximation. This corresponds 
to adding the contributions 
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k 4 \6(/) k \ 2 dk, 



(3.4) 



to Eqs. (3^) and ( p.3|j (see Refs. JlJJl^] for details). Note that this approach neglects mode-mode coupling and 
rescattering effects |14[ , which are important around the end of preheating, and the backreaction effect of metric 
perturbations. The Fourier transformed, linearized Einstein equations for field and metric perturbations in the 
longitudinal gauge are |34| 



(3.5) 
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d 2 U 



fc =2( 



(3.7) 



The analytic form of the solutions of the above equations are known in both the limits of fc — > and fc — > oo 
|Q. During inflation, metric perturbations exhibit adiabatic growth, <j>fe ~ cH/H 2 , after the first Hubble crossing 
(fc < aH). During reheating, the super- Hubble modes (fc <C aH) are nearly constant in the single field case. In 

contrast, the solutions for the small-scale modes (fc ^> aH) can be described by ~ (j> (cie lkv + C2e~' lkv ), which 
shows adiabatic damping during reheating. 

The system of scalar metric fluctuations $k and inflaton fluctuations S(f>k can also be described in terms of the 
Mukhanov-Sasaki variable Qk 0], defined by 



which satisfies the equation 



Qk + 3 H Qk + 



Qk 



H 



7<$>k 




(3i 



= 



(3.9) 



The study of this equation is particularly convenient, since, contrary to the equation of motion for $^ alone, 
nonsingular during the oscillations of the inflaton field. The gravitational potential is then related to Qk by 



it is 



iZL^! (En 
M 2 H \ rh 



(3.10) 



During inflation, modes of cosmological interest initially (at t = to) satisfy fc 3> aH, and their equation of motion is 
that of a free field in an expanding Universe. Quantization of this last quantity is thus straightforward. In the initial 
vacuum state only positive frequency waves are present (see for more details), so that we take as initial conditions 



Qk (h) 



1 



a (to) fcV2 



Q k (t ) = 



- i k 1 ' 2 

a 2 (to) 



(3.11) 



where ao is an arbitrary phase. 

As long as the modes remain much smaller than the Hubble scale they evolve as plane waves. When their physical 
lengths grow to of order the Hubble radius (fc ~ aH), the variation of the frequency of Qk becomes important. In 
fact, resonant amplification of some particular modes can occur during the coherent inflaton oscillations at reheating 
|l6f . This phenom enon is mostly appreciated when the inflaton is nongravitationally coupled to other fields (2(| (see 
subsection [II B 3 ). However, some resonant amplification occurs also in the single-field self coupled case with potential 
V = A0 4 /4 [p~8|, as we will discuss in the next subsection. 



B. Evolution of metric perturbations during reheating 

1. V = \m 2 (j> 2 

Let us first summarize the behavior of metric fluctuations during reheating in the single field massive inflationary 
scenario. During reheating, from the time-averaged relation (cp 2 } = (m 2 (j) 2 } one finds <j) ~ M p /(\/3tt mt) smmt . 
Neglecting in eq. (|1]) the H 2 and H terms, which decrease as ~ t~~ 2 during reheating, the rescaled variable Qk = 
a 3 ^ 2 Qk evolves according to a Mathieu equation with time-dependent coefficients 



dz 

where z = mt/2 , and 



— Qk + (A k - 2q cos 2z)Q k = 0, (3.12) 



W q= VT z - (3 - 13) 
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The small k modes (k < ma) lie in the resonance band around Ak = 1. However, the cosmic expansion makes q 
smaller than unity already after one inflaton oscillation, and the resonance is not efficient. As a consequence, modes 
in the long wavelength limit (fc — > 0) grow as Qk <x a 3 / 2 |4^-0], which makes Qk and nearly constant during 
reheating (for the mode k ~ 0.1m one has |fc 3 / 2 $fc| ~ 10~^7^ On the contrary the sub-Hubble modes (k > m) 
decrease due to adiabatic damping during reheating and are typically smaller than the modes which already crossed 
the Hubble scale by more than one order of magnitude. This behavior is shown for various wavelengths in Fig. |l|. 
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FIG. 1: The evolution of metric perturbations for k 



k/m 



0.1, 1, 5, 10 in the model V = m <f> /2 with inflaton mass m = 10 M p . 
We start integrating with an initial value, (j> = Q.5M P . The scale factor 
is normalized to unity at this time, and aH is always of order m in the 
time range plotted. Modes with k < m exhibit adiabatic growth, and 
then (when reheating starts, i.e., at mt m 1.5 in the plot) <J?fc becomes 
nearly constant. Sub-Hubble modes (k > m) instead decrease during 
reheating due to the adiabatic damping of the <j> term. 
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Let us next consider the single field massless chaotic inflation model, with potential V(4>) = A0 4 /4. As shown 
in (l^Q,^], the evolution of 4> is described (on dropping the term from the equation of motion, eq. ( |3.3| )) by 
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(3.14) 



where x = y/X tpo rj is the dimensionless conformal time, the suffix indicates the value of the quantities at the 
beginning of reheating and cn(x — xq, 1/v2) is the elliptic cosine function. 



Using the time-averaged relation 



\<p , one finds (H / H)' oc 



oc a 



and H cx a 



Hence, the Hubble 



terms in Eq. (3.9) rapidly become negligible. Using x and rescaling Qj. = aQk, Eq. (3.9) is reduced to 



dx 2 



Qk 



3cn^ 



V2 



Qk = 



(3.15) 



where we have defined k 2 = k 2 / (Ac^q) , and set xo = for simplicity. Eq. (3.15) is a Lame equation. From the 
analytical study of Ej| we deduce the presence of one resonance band in the interval 



-<*-< 



(3.16) 



G 



characterized by an exponential growth Qk oc e^ kX . The maximum value of the growth rate, /z max ~ 0.03598 , occurs 
at K 2 w 1.615. This growth dominates over the dilution from the expansion of the Universe even for the unrescaled 
variable Qk , and hence also for $^ . The growth of continues until the backreaction of inflaton fluctuations shuts 
off the resonance. The ma xima l value of the metric perturbation is found numerically to be \k 3 / 2 &k\ — 5 x 10~ 5 . 
Modes outside the interval ( 3.16 ) do not exhibit nonadiabatic growth, unless rescattering effects are taken into account 
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FIG. 2: The evolution of $ fc for k = 0.01, 1.25, 10 in the single field 
potential V = A<^> 4 /4 and A = 10~ 13 . The initial values of the numerical 
evolution are <f> = 0.1M p ,a = 1. We include the backreaction effect 
of produced particles via the Hartree approximation (see the text for 
details). Notice the enh ance ment of the mode K = 1.25 in the sub- 
Hubble resonance band (3.16). 



S.V= ±A<^ + ±sW 

We now discuss the situation in which a massless inflaton </> is nongravitationally coupled to another scalar field 
X, with conformally invariant potential V — A0 4 /4 + g 2 (f) 2 x 2 This model was studied in detail, in the absence of 
metric perturbations, in Refs. [ ff^]45| , showing the presence of an infinite number of strong resonance bands for x as 
a function of 



R = g 2 /X 



(3.17) 



The conformal invariance allows exact Floquet theory to be used to prove analytically that metric fluctuations can 
exhibit parametric amplification p6| . This makes the model a very useful testing ground for a number of issues 
p7j [30f including the possible formation of primordial black holes in some parameter space regions |2q| . Here we 
briefly summarize those of the above results which arc relevant for the calculations of the next section. 

As before, it is convenient to introduce a Mukhan ov-S asaki variable for the each field, Q k ^ = Sip k 1 ^ + (ip^ /if) 
i — <j),X which satisfy a simple generalization of eq. ( |3.9D p6| , p7| . 

Using eq. ( 3.14 ) and rescaling all variables as F = a F , the equation for Q k reads [^6| 
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where 
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(3.19) 



Here prime denotes a derivative with respect to ?y. In the early stages of preheating, as long as x fluctuations are small 
relative to 0, one can neglect the right hand side of eq. ( 3.18j ), which then reduces to the generalized Lame equation 
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studied in [ |44[|45| ] (in the following, we will apply the general analytical results for the Lame equation given in ]45fl). 
Of particular cosmological interest are the longwave modes on the scales proben by the current CMB experiments. It 
can be shown (45| that the resonance bands extends to very longwave modes (n — ► 0) only when the ratio R is in the 
range 



n (2n- 1) < R < n (2n+ 1) 



(3.20) 



where n is a positive integer. In each of these intervals, the growth rate [ik for Q£ can be analytically found as 



fJ-k 
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(3.21) 



where e = ^/2/Rk 2 and T ~ 7.416 is the period of inflaton oscillations. Then we have maximal growth when R is 
exactly 2n 2 and k = 0, with large characteristic exponent ^ max = (2/T)ln(\/2 + 1) — 0.2377. 

It is important to recognize that the exponential growth of just discussed does not necessarily translate into a 
parametric amplification of the gravitational potential 

47T 



$ fc + H<I> k = 



M 2 



In the two-field case, Eq. (3.5) is modified to 



(3.22) 



When x and Sxk are vanishingly small relative to and 8<pk, the evolution of $k is similar to the single field case. In 
fact, when R^$> 1 the homogeneous mode x as wen as the longwave Sxk modes (k — > 0) are exponentially suppressed 
during inflation |2^-|2^] due to the large effective mass <?0 relative to the Hubble parameter (see also Ref. fl46|- 148| ) . 
This can be explicitly seen by considering the equation of motion for Sxk , which for super- Hubble modes reads 

5'xk +3HS Xk +g 2 <f> 2 Sxk^0- 

(3/2-u) 



Hence, during inflation, the amplitudes of the long-waves modes evolve as |26| 



(3.23) 

\Sxk\ oc a-^/^-v) t w ith v = 

1/2 

Re [9/4— 3RM 2 / (27T0 2 )] . This confirms that the inflationary suppression of the long- waves modes becomes 
more important as R increases. In addition, considering backreaction in the Hartree approximation, one can show |2S(] 
that for R > 8 the growth of field perturbations shuts off the resonance before the super-Hubble metric perturbations 
(k — > 0) begin to be amplified. Hence, for large R , the evolution of the gravitational potential $fc on the scales 
probed by CMB experiments does not significantly differ from the "standard" one described in j34j unless rescattering 
i s tak en into account. For these modes, parametric amplification is expected to be effective in the first resonance band 
( |3.20 ) even when backreaction is considered |26|-|2l|. 
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FIG. 3: The evolution of metric perturbations <!>fc and field variances 
(Sx 2 ) and (8(t> 2 ) for « = 0.1, 3 in the potential V = A</> 4 /4 + g 2 tfx 2 /2 
with g 2 /\ = 2 and A = 1CP 13 . We start integrating with initial values 
cj> = O.IMp and x = 5 x 10 _6 A/ P , corresponding to x = 10" 3 M P at 
cj> — 4M P (i.e., 55 e- folds before the end of inflation). The final results 
are weakly dependent on the choice of the initial \ as long as 10 _6 M P < 
X < M p at <t> = 4M P ||,||. 
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In Fig. H we plot the evolution of <E>fc for the two modes n = 0.1 and k = 3 in the R = 2 case. We find that metric 
perturbations begin to grow after the 5\k fluctuations are sufficiently amplified. The growth of <J>fc ends around x = 90 
where the backreaction terminates the resonance, after which perturbations reach the plateau region found in Fig. ||. 

While only long-wave modes are necessary when comparing with CMB measurements, particle production is also 
induced by modes with larger momentum. The sub-Hubble modes (k > aH) at the beginning of reheating are 
not severely suppressed during inflation, and therefore the modes 6\k on those scales can be amplified to large values 
during preheating. However, from m x — g ((f)) we see that large values of R = g 2 /A results in a suppress ion o f 



the homogeneous mode \ during inflation. This also suppresses amplification of <£> from the rhs of equation (3.22) 
Indeed, for R ^> 1 the dominant effect for the source of metric perturbations is of second order in nonadiabatic 
(isocurvature) field perturbations 47 1. To estimate this growth, we consider the curvature perturbation on uniform- 
density hypersurfaces p9|,|50[] , denoted by ( . This quantity is related to the gravitational potential $ via f>0| 



C = K 



2p 



3(p + p) \aH 



(3.24) 



where 1Z = $ — (H/H) ($ + H $) is the comoving curvature perturbation, and p and p are the pressure and the energy 
density, respectively. The variation of £ occurs on large scales in the presence of a nonadiabatic contribution to the 
pressure perturbation, 5p na d — p(8p/p — Sp/p). One can estimate the second order contribution to C, denoted C(2)i 
as M 



3H 



2) 



5p na dHdt 



1 



2A^\ 
3iJ0/ 



g 2 (j) 2 8x 2 Hdt, 



(3.25) 



where we used the time averaged relation, 5\ — m^Sx — g <$> "X ■ We find that the second order term, g <jr5x , 
can induc e a va riation of while the first order term including the homogeneous \ is vanishingly small for R ^> 1. 
From eq. (3.25) we find 



V2ir 



2\(j> 3 \ 



) 



9 2 ^\ PsxlHdt, 



(3.26) 



where Ps x 2 k * s defined by = ^-(l^xll 2 )- Following Ref. j|7), the power spectrum of the 5\k fluctuation is 

estimated as 



Ps 



Xk 



2tt 



(3.27) 



Here kq ~ 0.15 is the mode corresponding to the horizon size at <p — 0.1Af p . From eq. ( 3.27 ) one finds 



^/|A V 2|2\ ^ t n x (|k1)n x (|k-k'|) 



H 



Hq\ 4 ( K \ 3 M(k,x) 



2tt 



K 



(3.28) 



where 



M(k,x) 




dn' / d6e 2 ^K' + ^-K>)* K > 2 s [ n e 
ii Jo 



dn'e^'^' 2 + 0(k 2 ). 



(3.29) 



Note that we expanded the term /j, k - k > around p K > in eq. ( |3.29 ). Due the k 3 factor in eq. ( 3.28| ), the spectrum P$ x i 
vanishes in the large scale limit (k — > 0), which implies that also second order effects in field perturbations do not 
induce large metric perturbations on the scales probed by CMB measurements pTj. In contrast, on sub-Hubble scales 
at the beginning of preheating (k ~ kq), Ps x 2 k IS nonvanishing. In this case parametric excitation of the x fluctuation 
can lead to a growth of ( and <&. 
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FIG. 4: The evolution of the gravitational potential and the curvature 
perturbation for k = 0.1 in the model of V = \(j> 4 /4 + g 2 (j> 2 x 2 /2 with 
g 2 /\ = 5000 and A = 10~ 13 . Note that we include the second order 
effect in field perturbations. We find that metric perturbations are 
enhanced on sub-Hubble scales. We also plot the evolution of the field 
variance, {8\ 2 ) and {5(f) 2 ), within the Hartree approximation. 



As an example, we plot in Fig. ^ the evolution of |fc 3//2 4>( 2 )| and |& 3 / 2 C(2)| for k = 0.1 in the case of g 2 /X = 5000. 
|fc 3 / 2 $(2) | is evaluated by making use of eqs. ( |3.24|) and ( |3.26|) . Note that the final field variances are smaller than in 
the g 2 /X — 2 case, due to the stronger backreaction effect (c.f. Fig. ||). The growth of nonadiabatic perturbations 
continues until backreaction ends the resonance. Although the final value of \k 3 / 2 Q( 2 ) | is somewhat smaller than that 
of \k 3 / 2 $ k \ in the g 2 /\ = 2 case, the fermion number density ( |2.16| ) can be larger due to the contribution of higher 
momentum modes, as we will see in the next section. 



IV. NUMERICAL RESULTS FOR THE PRODUCTION OF FERMIONS 



A. Quantum to classical transition 



Th e results of section [0] refer to particle production by a classical external source. On the contrary, the relation 
( |3. 10 ) of the previous section relates the modes of the decompositions 



Hv, x ) 



d 3 k 



V2 J (2tt) 



3/2 



[$ fc (v) e ik x a(fc)+^ (77) e- lk ' x at(fc)] 



(4.1) 



( r ?' x ) = -lx 



d 3 k 



y/2 J (2tt) 



3/2 



^ ( V ) e* x a(fc)+Q|« (77) J (k) 



(4.2) 



to the annihilation and creation operators a(k) and a'(k). As is well known, metric perturbations arise from quantum 
fluctuations during inflation and some of t hem e ventually undergo a transition to classicality. A consistent treatment 
of fermion production requires the integral ( |2 . 1 6| ) to be limited to those modes of $ that have grown enough to become 
classical. 

A simple criterion for classicality is given in Ref. pl| . Classicality is expected to be a good approximation when 
the following condition is satisfied: 



A 



W _ 



(i) 



» 1 



<t>,X 



(4.3) 
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The quantity is related to the uncertainty in the determination of (AQ^ AIlW) , where 11^ is the conjugate 



momentum of . In the case of pure vacuum fluctuations, the equality A)!' = 1 holds, in accordance with th e ru les 



(i) 



of quantum mechanics. Notice that this is the case for the initial states ( 3.11 ). A much bigger uncertainty (4.3) is 
associated to large fluctuations of classical nature. 



1 o° 



1 o° 



1 



1 o z 



1.0 t 



0.010 



1 



1 5 
t 



1 1 

r 






^ -q 




J k=a H ; 




V k=10a H \ 






, , , , i , , , , 


k=100a H i 

, , , , i , , , , i , , , , i , , , , ■ 



2 



2 5 



3 



FIG. 5: The quantum to classical transition as signaled by the evolu- 
tion of the quantity A^ during inflation and subsequent reheating stage 
in the single field model, V = A^ 4 /4. The mode k — aoHo leaves the 
Hubble scale at cj> = 1.5M P (i.e., at the initial time plotted). Af. grows 
significantly from unity after Hubble radius crossing. In this figure the 
end of inflation corresponds to t = \f\<j)ot ~ 10. 




FIG. 6: The spectra of A£ when the variance {8x 2 } reaches its maxi- 
mum value in the two-field model of preheating, V = A</> 4 /4+(/ 2 2 x 2 /2, 
for different values of R — g 2 /A. The small k modes (k < 0.5) are suf- 
ficiently enhanced by parametric resonance for R — 2n 2 , which makes 
the perturbations highly classical. In the case of R < 20 we find that 
A£ is very close to unity for the modes, k > 3. Therefore the cut-off 
k c = 3 is justified in evaluating the occupation number of fermions for 
R < 20. 
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In Fig. U the time evolution of is shown in the model V — A</> 4 /4, for some particular values of k. We 

notice that A^ begins to increase around the region where a mode leaves the Hubble scale during inflation. This 
is expected to occur quite generically, irrespective of the specific form of the inflaton potential. On the contrary, 
modes which are always inside the Hubble scale typically conserve their quantum nature ( A^ — 1 ) , unless some 
physical (model dependent) process drives them classical. Among these processes, nonperturbative particle creation 
during preheating can drive to a classical level the sub- Hubble modes in the resonance bands discussed in the previous 
section ^|H,|§|o|. 



A£ in the two-field model of preheating ( III B 3| ) for some values of R 



In Fig. H we show instead the spectrum 
(the largest time plotted corresponds to the moment at which Hartree approximation breaks down). The case R = 3 
reproduces the results of the one massless field case upon identification of x with the fluctuations of the field 4> ■ F° r 
this value, Fig. |^ shows the enhancement of the modes in the resonance band (3.16). In the two-field case, strong 
resonance occurs in the parameter ranges described by Eq. (|3~20|) . In the center of the resonance bands, R — 2n , 
small k modes (n < 0.5) are efficiently excited, while the modes with k > 3 are regarded as quantum fluctuations 
(A? ~ 1) as long as R is not much larger than unity (this evolution does not include rescattering effects which become 
important in the nonlinear regime). 

In the next section we numerically evaluate the formula for the occupation numbers reported in Sec. II. In the 
calculation we do not include those modes whi ch ar e not exc ited to a classical level (A[, 1 ' ~ 1). For modes which 
become classical, from the two decompositions (2.17) and (4.1) we get 



$(fc) 



V2 (2tt) 3/2 



(4.4) 



B. Particle production 



We can combine eqs. ( 2.16 ) and (4.4) to evaluate the number density of produced fermions in terms of the modes 
$fc computed in the previous section. Moreover, we assume the metric perturbations to be statistically isotropic and 
homogeneous, so that <E>k = $|k| ■ In this case, the final particle density is given by 



N2 
V 



^ 5 / d P0 / d PP 6 
15 (2 7r) or Jo Jo 



if 



dr)e ipor > <S> p (rj) 



(4.5) 



As discussed in the previous subsection, the integral over p = |p| must be restricted to the region where the per- 
turbations are classical. We have seen in section III that this typically occurs for p smaller than a cut-off value p c 



which approximately corresponds to modes which never crossed the Hubble scale (or, in the massless inflaton case, 
to the highest resonance band excited). In the numerical computation of N 2 we also integrated po up to p c . Higher 



frequencies are not expected to give a sizable contribution to eq. (4.5), since the rapidly oscillating phase e lpori averages 



to nearly zero the time integral. 

In the numerical calculations we take as the upper limit rjf of the time integral the moment in which backreaction 
effects start to shut off the resonance. In the single self-coupling case this occurs at Xf = VX(p(t co )r]f ~ 500, while in 
the two-field case the precise value of 77/ is a function of the ratio R . Note that this is just an approximate calculation, 
since the above integral is defined rigorously only for the asymptotic future, 77/ — >■ oo, with vanishing perturbations 
0>k —>■ . That is, our numerical calculation is approximately 



rty = a 3 N^ 



1 



15 (2tt) 



dpo / dpp 



if 



drje^Qp (77) 



(4.6) 



Before presenting the numerical results, we bri efly consider the occupation number of massive fermions produced 
by the homogeneous FLR W e xpansion, eq. ( [2.11 ). For sufficiently small masses, the final occupation numb er is still 



( appr oximately) giv en by fl4.q), since the Bogolyubov coefficient f3 is very close to zero in this limit (see eq. (2.8) and 
( 2.13 )). The term ( 2.11 ), although quadratic in f3 , is not suppressed by the small perturbations. Thus, by solving 
the equation 



N (tjv) = N 2 (tjv = 0) 



(4.7) 



we can estimate up to which mass the productions by the inhomogeneities is comparable with that from the 
homogeneous expansion. 



12 



The calculation of Nq has been performed numerically in 52 (see also references therein). For low fermionic masses, 
the gravitational production from the homogeneous expansion mainly occurs when the Hubble expansion rate equals 
to,/, . After this time, the comoving fermion number density is given by [ p2| 



horn 



Nq/V 



C a Tn~L ( 

TO^ 



3 a 



(4.8) 



In this expression, Hq denotes the value of the Hubble rate at the initial time io , where the scale factor a (to) is 
normalized to unity. The parameter a is instead the exponent appearing in the expansion law a (t) oc t a when 
H (t) — m^p (as is well known, the coherent inflaton oscillations give effective matter domination, a = 2/3, in the 



V 



m 



2^2 



/2 case, and effective radiation domination, a = 1/2 , for V = \<j) /A). 



the two values C2/3 — 3 x 10 and C1/2 — 10 have been numerically found 



For the first coefficient of eq. (|hq), 
3. 



1. V = ±m 2 2 



Let us fi rst present the numerical results for the massive inflaton case. The precise value deduced from the in- 
tegral (4.6) is sensitive to the cut-off p c . To discriminate which cut-off should be taken, we consider the quantity 
A^. We found that approaches unity around k = IOto, so we chose p c = 10m. In this case, for the comoving 
occupation number (normalizing a = 1 at <f> = 0.5 M p ) we have numerically found 



n,j, 



— „3 



a 6 N, k 



3 x 10 _14 to 3 



(just to give an example on how the final result is sensitive to p c , we found ~ 1 x 10 14 to 3 for p c 
n^~lx 10~ 13 to 3 for p x = 20to). 



(4.9) 
5to and 



Equating (h6) with (4^), the mass below which the production from inhomogeneities dominates over the one 
from the homogeneous expansion is estimated as 



TO^ 



n 



C2/3H0 



10 2 GeV. 



(4.10) 



Of more physical interest is the fermionic abundance deduced from (4.9). By assuming an instantaneous inflaton 
decay into a thermal bath characterized by the reheat temperature T r h, the particle density divided by the entropy 
density, s, is given by 



(T rh ) = ^ (T rh ) 
s 



10- 



10 9 GeV 



(4.11) 



In the case of other light gravitational relics, such as gravitini or modulini, the primordial abundance is severely 
constrained by the successful predictions of Big Bang nucleosynthesis. For masses of order the electroweak breaking 
scale, the limit is about Y < 10~ 13 - 10~ 12 p3|. In the present case, eq. (4.11) gives a much smaller abundance. 



2. V = ±A0 4 



In the quartic case, both super-Hubble modes and the modes in the inte rval ( |3.16 ) contribute to the particle 
production, eq. ( |4.6| ). From simple estimates (for example combining eq. (|4.6|) with the values for the perturb ation s 
shown in Fig. |^), one can see that the main contribution to particle production is given by the resonance band (3.16). 

In calculating the total number densities of produced fermions due to inhomogeneities, we take the cut-off value in 
momentum space at p c — 3, over which A^ approaches unity. Then the final particle density reached at the end of 
preheating is numerically found to be 



1 x 10" 



(v^o) 3 



(4.12) 



which is larger than in the massive inflaton case. 

Proceeding as in the previous subsection we find that the production from metric perturbations dominates over 
that from the homogeneous expansion as long as the fermions have masses smaller than 



13 



1 



Ho \C\/2 



2/3 



10 5 GeV 



(4.13) 



The final abundance is also greater than in the previous case (thus confirming that the integral (4.6) is dominated by 
the modes in the resonance band), 



r^(T rh )~io- 22 

Nevertheless it is still too small to have any cosmological effect. 



(4.14) 



3. V = + \g 2 4> 2 x 2 



Since this potential leads to an infinite number of resonance bands it is the case of most interest. In Fig. ^| our 
numerical results for fermionic p rodu ction are presented as a function of R = g 2 /\ , for R < 10 . This is obtained by 
solving the linearized equation (3.22) using the Hartree approximation for field pert urba tions. By c omparison with 
the single field case (g = 0, here reproduced by the choice R = 3) the two equations (4.13) and (4.14) are generalized 
to 



rrijp 



5 x 



10 13 n 2 / 3 GeV 



(T rh ) ~ 10 



9 n-^ 



(4.15) 



where = n^,/(yA </>o) 3 • 

The qualitative behavior of Fig. is readily understood from the structure of resonance (3.2C). In the range of R 
plotted, we have maximal production when R is in the center of the first resonance band, in which case n$ is found 
to be ~ 5 x 10~ 10 , which is about 5000 times larger than in the single-field case. In the second resonance band, 
particle production is slightly smaller than in the R = 2 case. 




g 2 /X 



1 o 



FIG. 7: The comoving number density, n$ = n^j (\/A<^o) 3 , vs the ratio 
g 2 /\, for ff 2 /A < 10. Note that this is based on the linear calculation 
using eq. (3.22), which does not include the second order effect in field 
perturbations. We take \ = 10 _3 M P at <j> = 4M P . We see t hat p article 
creation is enhanced in the super-Hubble resonance bands (3.2C). 



As R increases, the inflationary suppression of x becomes more and more relevant, and particle production from first 
order perturbations is consequently reduced. Nevertheless, the second order effect of field perturbations described 
in subsection III B 3 can lead to the excitation of metric perturbations on sub-Hubble scales. Just for indicative 
purposes, we report a couple of our numerical results in this range. For large R, the resonance bands cover only 
modes up to |45[| 




(4.16) 
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In the center of the 10-th instability band, R = 200, one finds resonant amplification up to k < 1.78, as can be 
appreciated in Fig. ^. By implementing the analysis reported in subsection III B 3 we find n^p ~ 6 x I P -11 . As R 
increases, higher momentum modes contribute to the growth of metric perturbations. For example, eq. (4.16) gives 
k < 4 in the R = 5000 case. As a consequence, this enhances fermionic production. The comoving occupation number 
in this case is found to be ~ 1 x 10~ 7 , Y$ (T r h) — 10~ 16 . 

It is however important to mention that for 7? 3> 1 the whole analysis becomes very delicate. For large coupling g , 
the backreaction effect of produced particles ends resonant amplification of field perturbations earlier. In addition, 
it was found in Ref. (li) that the final field variances get smaller if rescattering of the 5\ fluctuations is taken into 
account in the rigid FLRW spacetime ($ = 0). On the contrary, the effect of rescattering can lead to the excitation 
of inflaton fluctuations through the amplification of 5\ fluctuations. Since (contrary to the field x) the homogeneous 
inflaton component is not suppr essed during inflation, the gravitational potential can then acquire a potentially large 
additional source (see eq. ( |3.22| )). Finally, the anisotropic stress is expected to be important in the nonlinear regime, 
in which case the relation $ = ^ no longer holds |3(]] . The study of these effects is certainly worth separate detailed 
investigation. 



V. CONCLUSIONS 



In the present work we have discussed gravitational creation of fermions by inhomogeneous scalar perturbations 
of a Friedmann-Lemaitre-Robertson- Walker (FLRW) Universe. Massless fermions are conformally coupled to this 
background, and thus are not generated by the homogeneous gravitational field. As a consequence, gravitational 
production of light particles is solely induced by metric perturbations which break the conformal flatness of the 
background. Although the small size of perturbations on very large scales indicated by CMB experiments seems to 
suggest a small production, a precise computation is still worth while for light gravitational relics. 

This is particularly true in scenarios in which metric perturbations undergo parametric amplification during pre- 
heating at the end of inflation. To examine this issue more quantitatively, we have studied some particularly simple 
models of chaotic inflation. The simplest example is a single massive inflaton field (V = m 2 4> 2 /2), where parametric 
resonance is actually absent. In this model, metric perturbations are nearly constant during reheating for modes which 
left the Hubble scale during inflation, whereas the modes deep inside the Hubble scale exhibit adiabatic damping. As 
may be expected, particle production is negligibly small in this situation. 

We then studied the massless self-coupled inflaton field (V = A0 4 /4) case, in which parametric resonance occurs for 
modes in a small band in momentum space near the Hubble scale. As a consequence, particle production is enhanced. 
Moreover, for a massless inflaton the background energy density decreases more quickly with respect to the massive 
case, and this results in a larger abundance at reheating of the particles generated by the inhomogeneities. Despite 
these two effects, the final result for the production is however still well below the limit from nucleosynthesis. 

More efficient production is expected when more fields are present, such as the model characterized by the potential 
V = A</> 4 /4 + g 2 2 x 2 /2. In this case, the crucial parameter which determines the resonance bands is R = g 2 /X- 
In particular, the longwavelength modes of the second field x are parametrically amplified whenever R ~ 2n 2 with 
integer n . At linear order in field perturbations, these modes are coupled to metric perturbations through a term 
proportional to x$Xk ■ F° r relatively small R , metric perturbations are amplified by this term, and particle production 
is strongly enhanced. 

For R > 10 the large x effective mass (m x — g{4>) <x R 1 / 2 ) causes the \ ne ld to be exponentialy suppressed during 
inflation. As a consequence, for large R metric perturbations mainly increase due to second order field fluctuations Ji7| . 
Particle production may be expected to become more significant at very large R . This is indeed suggested by our 
numerical results, although in the cases considered the final abundance is still below the nucleosynthesis bounds. 

Our numerical simulations include the backreaction of the fluctuations on the background evolution in the Hartree 
approximation. Rescattering effects - mode-mode coupling between fields - are instead neglected. This may consid- 
erably affect the final abundances, particularly in the large R limit where the mode coupling becomes particularly 
important. Strong rescattering effects imply an earlier end to preheating, which may decrease the final values of field 
and metric perturbations. On the other hand, rescattering leads to amplification of modes outside the resonance 
bands. In particular, inflaton field fluctuations can be significantly amplified. Since (contrary to the field x) the 
homogeneous inflaton component is not suppressed during inflation, metric fluctuations may be strongly enhanced 
by the source term proportional to 4>6<pk (analogous to the x$Xk term considered above). This would be reflected in 
more particle creation. 

Stronger production from inhomogeneities may also occur in different contexts. In the multi-field inflationary 
scenario, the gravitational potential can exhibit nonadiabatic growth during inflation |54|,|55|], especially when the 



second field x nas a negative nonminimal coupling to the Ricci scalar [p6 57 . In addition, it was shown that the 



15 



growth of metric perturbations can be present even in the single field case in the context of non-slow roll inflation 
with potential gap |5§f| , and during graceful exit in string cosmology fo"^] . It may be interesting to study gravitational 
particle production in these contexts along the lines presented in this work. 
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